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In this paper we derive general relations for the band- structure of an array of quantum dots and compute its 
transport properties when connected to two perfect leads. The exact lattice Green's functions for the perfect 
array and with an attached adatom are derived. The expressions for the linear conductance for the perfect array 
as well as for the array with a defect are presented. The calculations are illustrated for a dot made of three atoms. 
The results derived here are also the starting point to include the effect of electron-electron and electron-phonon 
interactions on the transport properties of quantum dot arrays. Different derivations of the exact lattice Green's 
functions are discussed. 
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I. INTRODUCTION 

Quantum dots, quantum wires, and molecular structures 
are among the most studied low-dimensional condensed mat- 
ter systems due to their importance to nanoelectronics, ' and 
recently to biology.^ In the field of molecular electronics, 
molecules are used to control the current flow when they 
are assembled in between metal contacts.-' Among the sev- 
eral physical properties exhibited by low-dimensional elec- 
tron systems, interesting ones can be found in the transport 
properties of these systems, such as conductance quantization 
and conductance oscillations,"* the latter effect depending on 
the number of atoms in the wire. Some of these effects have 
been found in the related field of carbon nanotubes, where ex- 
periments have shown that the conductance of a single wall 
carbon nanotube is quantized*" and shows Fabry-Perot inter- 
ference patterns, a signature of coherent transport in the car- 
bon wire. 

The above systems all fall under the study of low- 
dimensional physics and are best described starting from the 
tight-binding approximation. When we address the physical 
properties of these low-dimensional systems there are two dif- 
ferent questions to address. One is concerned with the elec- 
tronic spectrum,- and the other with the transport properties.^ 
The first of these properties is of importance for the optical re- 
sponse of the system and the second for its use in nanoscopic 
devices. 

The need of an efficient method of computing both the spec- 
trum and the conductance of these systems is therefore obvi- 
ous. Both these problems can be traced back to the calcu- 
lation of the lattice Green's function of the system. In the 
case of the spectrum, the Green's function should be com- 
puted for the isolated system, whereas when considering the 
transport properties, the Green's function must be computed 
taking into account the effect of the coupling between the 
system and the metallic leads. The lattice Green's func- 
tion method was used to describe the appearance of surface 
modes (Tamm states) in a finite one-dimensional chain and 
its interaction with a non-linear impurity.^ The generalization 



to a semi-infinite square lattice was also discussed.^ Several 
different methods of computing transport properties of one- 
dimensional quantum wires, using lattice Green's functions, 
are available in the literature. A tutorial overview on some of 
the used methods was recently written by Ryndyk ef oL . '° Us- 
ing the Keldysh method, the coherent transport through a one- 
dimensional lattice was studied by Zeng etal.J^ Li et at. stud- 
ied the transport through a quantum dot ring with four sites. — 
The inclusion of time dependent potentials on the transport 
properties of one-dimensional chains was done by AiTachea. — 
The inclusion of electron-electron interactions in the transport 
properties of a small system was discussed by Oguri'"* and of a 
quantum wire was done by Karrach et al.,^^ using a functional 
renormalization group method. The extension to a quasi-one- 
dimensional Kagome wire, where the feature of a multiband 
system is present, was considered by Ishii and Nakayama. - 
The interesting situation where the metallic wire is connected 
to a Heisenberg chain was studied by Reininghaus et at. , - The 
above results are just a very small subset of the existing rep- 
resentative literature on quantum transport, where the concept 
of lattice Green's function plays a central role. 

Historically, the first approach to electronic transport in a 
one-dimensional finite system was done in a series of elegant 
papers published by Caroli et at. In these works the au- 
thors addressed the question of how defects affect the charge 
transport of the quantum wire. Indeed, the question of how 
localized defects change the otherwise perfect transport prop- 
erties of the system has been addressed by several authors in 
the framework of tight-binding systems i^''^^'^^ 

Guinea and Verges^' used a Green's function method to 
study the local density of states and the localization length 
of a one-dimensional chain coupled to small pieces of a poly- 
mer. They showed that at the band center there is a complete 
suppression of the transmission coefficient due to a local an- 
tiresonance. Sautet and Joachim^^ studied the effect of a sin- 
gle impurity on the transport properties of a one-dimensional 
chain. The impurity was assumed to change both the on-site 
energy and the hopping to the next-neighbor atoms. Mizes and 
Conwell^^ considered the effect of a single impurity in two 
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coupled-chains, showing that a change on the on-site energy 
has a more pronounced effect in reducing the transmission in 
the one-dimensional chain than it has in this system. Finally, 
Peres and Sols studied analytically the effect of a localized 
defect in the transport properties of polyacene (a multi-band 
system), putting in evidence a parity effect, that was also used 
by Akhmerov et al?^ to formulate a theory of the valley-valve 
effect in graphene nanoribbons .^^-^^-^^ Also the study of va- 
cancies in transport of quasi-two-dimensional systems is an 
important area of research. The study of a linear array of 
quantum dots, represented by a single site with an s— orbital 
was carried out by Teng et al.}^^ 

The book by Ferry and Goodnick has a very good introduc- 
tory section to lattice Green's functions.^' Also the recursive 
Green's function method is there presented. The method we 
present in this paper has strong similarities to the recursive 
Green's function method;^ The important difference to stress 
here is that the recursive Green's function method is imple- 
mented as a numerical method, whereas our approach does 
solve the same type of problems in an analytical way. The link 
between Green's functions and transport properties of nanos- 
tructures and mesoscopic systems is well covered in the book 
by Datta.^^ 

In this work we give a detailed account of a method, based 
on the solution of the Dyson equation, to compute the lat- 
tice Green's function of an array of quantum dots, where the 
dot is represented by an arbitrary number of sites with arbi- 
trary values of the site energies and of the hopping parame- 
ters. The method is developed within the approximation that 
there is only one hopping channel between the dots in the ar- 
ray. This constrain is used to keep the level of the formalism at 
its minimum. The formalism is easily generalized to include 
the effect of defects (both on-site and adatoms defects) and to 
describe the transport properties of both the clean and the per- 
turbed system. The method exploits the fact that the Dyson 
equation for the lattice Green's function can be solved exactly 
for bilinear problems as long as the hoppings are not of arbi- 
trary long range. Also the developed formalism can be used to 
describe surface states and non-linear impurities effects^ in a 
finite array of quantum dots, but we will not pursue these two 
aspects in this paper. Although developed within a quasi-one- 
dimensional perspective, we will show in a forthcoming publi- 
cation how the method can be generalized to two-dimensional 
ribbons. 



11. DETERMINATION OF THE LATTICE GREEN'S 
FUNCTION 

In this section we develop a general method for determining 
the Green's function of an array of quantum dots. From it both 
the electronic spectrum and the rule for momentum quantiza- 
tion are obtained. As a warming up we first revisit the solution 
of the finite chain problem. 

The traditional approach^^ "''*'^^ to determine the Green's 
function in real space requires the previous solution of the 
Schrodinger equation, with the corresponding determination 
of its eigenvalues and eigenvectors. After this is done, the 



Green's function is computed using 

G(r,r',.) = 5:^^«M::i (1) 

where ^/'„(r) and A„ are, respectively, the eigenstates and 
eigenvalues of the eigenproblem Hipn{r) — AnV'nlr), where 
H stands for the Hamiltonian of the problem and the summa- 
tion over n means a discrete sum (for discrete eigenvalues), or 
an integration (for continuous eigenvalues), or both, over the 
quantum numbers of the problem. 

The evaluation of the sum in Eq. ([U may be a very hard 
task, depending on the mathematical complexity of both the 
wave functions and the eigenvalues. Even for relatively simple 
cases the evaluation of the summation is far from obvious (see 
appendices of Ref. [36l) . 

An alternative approach, used for lattice systems, starts 
from the definition of the resolvent operator 

G = , (2) 

z — H 

and computes the matrix elements of the resolvent by eval- 
uating directly a number of determinants associated with the 
matrix (z — H).^^'^'' This method has the obvious drawback of 
being limited by the possibility of computing analytically the 
necessary determinants.i^^^ 

We present in what follows a method that overcomes the 
technical difficulties above mentioned. 

A. The single chain case 

In order to understand how the method works, we revisit the 
problem of determining the Green's function of a finite one- 
dimensional chain of atoms, with a single orbital per atom. 
This simple example will help us to fix the notation and state 
the general arguments about the solution of this type of prob- 
lems. Let us assume that the system has N atoms, with the 
motion of the electrons described by the tight-binding Hamil- 
tonian 

H = Ho + V, (3) 

with 

N 

Ho = eoJ2\^){i\, (4) 

and 

Af-l 

v = -tJ2i\^){^ + M + \i + ^){^)- (5) 

i=l 

Clearly, Eq. (|4|i represents the on-site energy of the electrons 
in the atoms and Eq. (|5]) represents the hopping of the elec- 
trons between neighboring atoms. This may seem as an im- 
portant restriction, but it can in fact be relaxed and the ap- 
proach extended to more general hopping processes.'*" Alter- 
natively the Green's function of a more complex Hamilto- 
nian may be generated using the extension theory for lattice 
Green's functionsiii 
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Let us now introduce two different resolvent operators, the 
free resolvent 



1 



(6) 



and the full resolvent G, given by Eq. (|2]) with H given by Eq. 
(O. The strategy is to determine G by solving exactly Dyson's 
equation,^-* considering the hopping term V as perturbation. 
In terms of the resolvents and of V, the Dyson equation takes 
the form 



G ^G° + G^VG . 



(7) 



Forming matrix elements with the basis vectors and using 1 



N 



^ one obtains 
1 



GO 



^\G\J)=S.,,-t{{^-l\G\]) + {^ + l\G\J))., 



(8) 



with G" = (z — eo)^^ . Apart from the term 5i,j, this is the 
equation for the wavefunction of a particle of coordinate i with 
the tight binding Hamiltonian (|3]l, and eigenvalue G^{z). It is 
obvious that the difference of two solutions of this equation 
will be a solution of the equation without the 5^ j term. So 
G 



{i\G\i) can be determined by adding a general solu- 
tion of the homogeneous equation to one particular solution 
of the full equation. The latter can then be determined by the 
boundary conditions. 

The solutions of the homogeneous equation are superposi- 
tions of plane waves, {n\G\m) = Ame^^"' , where Am is an 
arbitrary function of m, and 9 is defined by 



1 

G" 



(9) 



the usual dispersion relation for a ID tight-binding problem 
with nearest neighbor hopping. To find one particular solution 
of the full equation, we use the fact that it has to satisfy the 
homogeneous equation for i < j and i > j and therefore be 
a linear combination of plane waves of wavevectors ±0 — the 
solutions of Eq.® — or, equivalently, of snv{id) and cos(i6'): 



A> cos{ie) + sin(i6l), for i < j , (10a) 
: G^ cos(i6i) + sin(i6'), for i> j . (10b) 



In the linear system of equations obtained by fixing j in Eq. 
® , 



1 <G0 0" 




\ (l|G|j) 1 







tG° 1 tG° •••0 




(2|G|j) 







tG° 1 iG" ••• 














{J\G\3) 




GO 


tG° I 




. {N\G\3) . 








(11) 

all but two equations are automatically satisfied by the fact 
that Gfj and G?^ solve the homogeneous equation, leaving 



only two conditions mixing G^ with Gf^ : 
1 



n< 

QO jj 



1 - t {Gf_^^j + Gf_^j) 



1 

GO 



(12a) 
(12b) 



These are easily shown to be equivalent (using the fact that 
the Green's function is also a solution of the homogeneous 
equations) to 



G< 



Gf] 
1 

t ' 



(13a) 
(13b) 



The Equation ( |13at corresponds to the continuity of the 
Green's function, whereas Eq. ( I13bl l corresponds to the dis- 
continuity of the derivative of the Green's function in the the- 
ory of second order differential equations.^"* 

Inserting Eqs.lfTOll in Eqs.(fT3]l , one can obtain a rather sim- 
ple solution (valid in both domains i < j and i > j) as 



Gij — 



1 



2t sin( 



:sin i9\i-j\) 



(14) 



The general solution is obtained by adding an arbitrary solu- 
tion of the homogeneous equation. 



1 



2t sine 



-J\ 



Gij = A^ cos{9i) + sm{di) 

/./, >i I I I 1/ 

(15) 

The free coefficients, A^ and B^ , are determined by boundary 
conditions. For a finite chain we must enforce. 



(0|Gb-)-(7V + l|G|j) = 
(i|G|0) = (i\G\N + l) = 0. 

It is straightforward to show that these lead to 

1 cos{N9 + 9) 



Gij (z) 



2t sm{N9 + 9) sin 6* 



[cos(i6' — jL 



1 

2t 



s'm{i9 + jO) sm\i9 — jl 



sin 9 



(16a) 
(16b) 



cos{i9 + j9)] 
(17) 



Noticing that the Chebyshev polynomials obey the finite dif- 
ferences equation 



fn-l{x) + fn+lix) = 2xf„{x) 

and have the representation^ 

Tn{cos9) 

Un{cOs9) 



cos{n9) , 
sm{n9 + 9) 



(18) 

(19a) 
(19b) 



we see that our solution (fTTT i is the same obtained for a finite 
one-dimensional harmonic lattice,"*^ as it should be. It is worth 
noticing that the rule for momentum quantization (that is 9) is 
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obtained from the poles of the Green's function ( fTTI l. leading 
to 



sin(iV + l)^ ^ 0:6i 



ni 



7V + 1 ' 



(20) 



with£ = 0,l,2,...,iV-l. 

Rewriting G° = (z - eo)"^ as G° = {E - €a + iO+)-\ 
where E stands for the energy of the electron, the energy de- 
pendence of the Green's function is obtained;^^ this allows the 
calculation of the local density of states from 



puiE) = --'iGu{E + iO+). 

TT 



(21) 



It is clear that pu (E) is site dependent for a finite chain. It 
is an elementary calculation to find an explicit expression for 
Pii(E), by writing both cos(2i6') and sin(2i6') in terms of sin 6 
and cos 9 and using Eq. (|9|l. We conclude by stressing that the 
matrix elements of the resolvent were obtained without the 
need of evaluating the sum over the eigenstates, as in Eq. ([T]i, 
which constitutes the major advantages of this method over 
most of the existing ones. 

As a final comment in this section, we remark that the final 
solution for Gy is quite symmetrical in both coordinates, even 
though the method we followed treats the two coordinates on 
a rather different footing. In Appendixes |A] and |B] we give 
two other derivations of the Green's function of a finite chain 
which treats both coordinates equally from the beginning. 



B. The general case 

We now consider the case where N quantum dots, repre- 
sented by hexagons in Fig. [T] are coupled together by an 
hopping parameter t. Although having just a single hopping 
channel between the dots may seem rather restrictive, we used 
it nevertheless to illustrate the method. Also, if we choose to 
position the dot in such a way that is has oriented edges such 
that a particular site of the dot is closer to the next dot than any 
other point, as in the case of Fig. [T] the used approximation 
is somewhat justified. The generalization of this possibility to 
multi-hopping hopping processes does not changes the gen- 
eral idea, but adds some complexity to the final solution. In 
addition, we assume that the dot is also described by a lat- 
tice model. Although this an apparent restriction, it can also 
be relaxed. We now address the question of determining the 
Green's function of the quantum dot array. 
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N 



Figure 1 : An array of A'^ quantum dots coupled together by an hop- 
ping integral t. The dot is represented by the hexagon. The points 
I and r are those that are coupled between dots due to t, and the 
choice of the letters stand for most left and most right points in the 
dot, respectively. 



The Hamiltonian for the problem defined by Fig. [T|is writ- 
ten as 



Ht = H° + V , 



were 



N 



H 



D 



(22) 



(23) 



with Hp the Hamiltonian of the i dot (which is not necessary 
to specify at this point), and 



V = -tJ2i\hr){i + l,l\ + \i + l,l){i,r\) 



(24) 



As before we define G'^ — {z — H^)~^, and note that \i, f3) 
are the basis states with i ~ 1,2, ... ,N and /3 labeling the 
sites in the quantum dot; we will need two different /3's only, 
f3 = r,l. 

The matrix elements of the Dyson's equation (O, given the 
Hamiltonian (l22b . reads 

{t,a\G\j,f3)^S.,,G% - t{G°^i{z-l,r\G\j,P) 

+ Gl.{i + l,l\G\j,P)),(25) 

with G'^p = {i,a\G'^\i, (3). It is easy to see, by direct replace- 
ment, that the homogeneous Dyson's equation (that is when 
i ^ j) is solved by an Ansatz of the form G^™ = A™^e*"^ 
when 6 is chosen such that 



1 + 2tGl cose + t'iGlr - rG'iiG';.^ 



(26) 



an expression that gives the band structure once the quan- 
tized values of 9 have been determined. We see from Eq. 
( |25l ) that {i, a\G\j, j3) is only coupled to (z — 1, r|G|j, j3) and 
(« + 1, l\G\i, (3). Therefore we will solve the Dyson's equa- 
tion dZSl) for the particular case of a, /? = I, r. To start with 
we make the Ansatz (a linear combination of terms of the form 

Am AnO\ 



Glp = A<pCO<^^)+B<pSm{i9),fov i<j, (27a) 
G^f = Alp cos{i9) + B>p smii9) , for i>j, (27b) 

where G^^'"* = {i, a\G\j, f3), and the multiplicative coeffi- 
cients of the trigonometric functions depending on j. The 
finiteness of the chain is imposed by the conditions 

(0,r|G|j, /?) = {N + l,l\G\j,P) ^0, (28) 
and time reversal symmetry implies that 

{i,a\G\j,P)^{j,P\G\i,a). (29) 

The solution of the non-homogeneous Dyson's equations (that 
is, when i — j) leads to a linear system of equations of the 
form 



Mv = b, 



(30) 
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with the transpose vectors given by = [Af'i , Af"^, B^^] ef a/. 3 The Hamiltonian of the dot we are considering reads 
b"^ ~ [G^, G^,^, 0] . The matrix M is easily constructed from 
the non-homogeneous Dyson's equations and is given in the 
AppendixICl The solution of the linear system is best accom- 
phshed by Cramer's rule, leading to 



A> 



qO 

" '■'■■sin(7V6' + 6l) sin(j6l) 



Pl 

rr " 

Pl 



Pn-j 



(31a) 
(31b) 
(31c) 



with Pl = tGl sin OPl and = tG^^ sm{xe~6)+sm{xe). 

Combining the homogeneous and non-homogeneous equa- 
tions we can derive the following results 



^<,m— Im ^>.m — Im 



r 



(32a) 



which combined with time reversal symmetry and one of the 
non-homogeneous equations leads to the linear system 



Vu = q, 



(33) 



withn^ = [A>,B>,B<] 



i^,.;, i^j^j, and V and q given in Appendix0 
The solution of the linear system gives 



A 



{tPL)-'PLPm-i , (34a) 

-{tPLy^[cOs{Le) + tG^ COS{L0 - 0)]PrrMh) 
P^r = {tPLY^[cOSe + tGl.COs{2e)]PL^ra- (34c) 



P^l 



-^A a)(i, a + 1| + |i, a + a|) 

Q = l 

3 

+eo^ |i,a)(i,a| , (36) 

a=l 

with the boundary condition (i, 4| — (i, 1| . What is now nec- 
essary is to compute the Green's function for this system. It 
just happens that the Green's function for this system is given 
by the Green's function of a finite chain (three sites) with pe- 
riodic boundary conditions. This can be obtained from the 
procedure of Sec. Ill Al by replacing the boundary conditions 
(I16ab bv 



which leads to 

Gil — 



1 



1\G\])^{N + 1\G\3). 
sin(iV6') 



(37) 



2t A cos(A^6') sin 9 - sin 9 
^ -[/|,_j|_i(cos6l) . 



T,.,{cos9) 



2tA 



(38) 



For the Hamiltonian Qoj we have = 3 and the Green's 
functions computed from d38T l are given by 



rr 



G 



Li,., 



II — "-^11 
^12 



1 



1 + 2A 



2tA -1 + A + 2A2 ' 

1 1 
2t^ 1 - A - 2A2 ' 



(39a) 
(39b) 



As we saw before the poles of the Green's function gives 
the rule of momentum quantization. In this case the poles 
correspond to the zeros of Pl, which leads to an equation for 
the quantization of 9 that reads 



tG'^^ sm{N9 -9)+ sm{N9) = . 



(35) 



Contrary to the case of the single finite chain, Eq. (|35] | de- 
pends on the energy, and therefore it has to be solved together 
with Eq. (ESJ. The knowledge of A>, A>, A>, B>, B<, and 
B^^ is all that is necessary to determine all the Green's func- 
tions for this problem. The full form of the Green's functions 
is given in Appendix iDl 



C. Explicit results for a particular example 

Let us now consider a specific example and work out the 
energy spectrum and the momentum quantization. For the 
Hamiltonian of the quantum dot we consider the case where 
the dot is made of three sites very close together, with a single 
local orbital^** associated to each site. The sites are coupled 
together by a hopping matrix element t^- The simpler case of 
representing the dot by a single site was considered by Teng 



with X = {E - eo)/(2tA). We should note that G?i can be 
written as 



1 



1 1 



2tA V3 A - 1/2 3 A + 1 



(40) 



which means that the eigenvalue A = —1/2 is bidegenerate, 
and that is the fundamental reason why the denominator of the 
Green's function is not a cubic polynomial. 

We shall now consider the physical relevant case where 
t A ^ t. Within this approximation the solutions of Eq. 
are given by (A ^ 1) 



-1 


t 

' 3tZ 


COS 9 , 


1 


t cos t 


? t 


2-*- 


6tA 


6tA 


1 


t COS t 


1 t 


2^ 


6tA 


^ 6tA 



(41a) 



A2 

A3 



The values of 9 are obtained from the solution of Eq. ( |35]) . 
which requires the knowledge of G52, which in the approxi- 
mation of Eq. (l4n i are given by 



6 



G?2(Ai) 



r 



2 cose 



r 



cos 6* + V3 + cos2 e 



COS 6* + V3 + cos2 61 
The Eq. ([35) is trivially solved for the case A = Ai, giving 



(42a) 
(42b) 
(42c) 



N +1' 



(43) 



with £ = 0, ±1, ±2, . . . , TV. The other values of 9 for A2 (top 
sign in Eq. (l44l i) and A3 (bottom sign in Eq. (l44l i') are obtained 
as solutions of 



i{Ne) sin 61 ± sin(7V6l) V3 + COS^ t) , 



which for large N reduces to 



(44) 



(45) 



with t ~ ±1, ±2, . . . ,N. This result can be appreciate 
graphically, by plotting both sides of Eq. ( |44| ) on the same 
graph. Naturally, when N 00, 9 becomes a continuous 
variable in the interval 9 e [— 7r,7r[. In Figure |2] we plot 
the eigenvalues A^, given by Eq. ( |4TI ). using t/^A = 0.1. 
The three sites composing the dot originate three energy mini- 
bands (since i a ^ t) in the dot array. 





Figure 2: Eigenvalues Ai, when i — 1 (left) and i = 2, 3 (right) 
as function of 9/tt, for large values of A'^. We have used the ratio 

t/tA=0.l. 

It should be now clear that this method allowed us to deter- 
mine the energy spectrum, the quantization rule for 9, in the 
case of a finite N, and the Green's functions for the array of 
quantum dots, with essentially the same effort it would take 
us to solve the Schrodinger equation for this problem. 



ni. THE TRANSPORT PROPERTIES 

We next want to work out the transport properties of the fi- 
nite system described in Sec. Ill CI We first derive the general 



results and latter use them to study the system introduced in 
Sec. Ill CI We will consider that our system is connected to 
two semi-infinite perfect leads. The leads to which we will 
connect the system will be described in similar terms to those 
used in the Newns modelSM'' and the connection between the 
Green's function and the transmission across the system is de- 
scribed using a formalism similar to that introduced by Mujica 
et flZ.,'*^''**''"^^ which, it turns out, it is similar to the approach de- 
veloped by Fisher and Lee.^" The general relation between the 
transfer matrix approach and the Green's function method is 
described in Ref. 1.44.1 . in the context of continuous models. 



A. The general formalism 

Next we present the general formalism including some par- 
ticular important results, associated with the concept of sur- 
face Green's function. We represent the Hamiltonian of the 
perfect leads by one-dimensional semi-infinite tight-binding 
models reading 



Hl = ~Pl {\^-'^,L){i,L\ + \t,L){z-l,L\X^6a) 

i — — OD 

00 

Hr - ~I3r {\i + l,R){i,R\ + \i,R){i + l,m^) 

i=N+l 

where /?l,_r are the hopping parameters of the left (L) and 
right (i?) leads. Although we are representing the leads by 
a one-dimensional model, this is of no consequences in the 
characterization of the transport properties of the dots, being 
only essential that Pl.b. are such that metal bands in the leads 
have very large band-width. Since the effective band width of 
our dot structure is proportional to t, the only condition is that 
Pl,r ^ t. The coupling between the leads and the dots is 
made by the Hamiltonians 

Vl = tL{\0,L)(l,l\ + \l,l){0,L\), (47a) 
Vr = tRi\N + l,R){N,r\ + \N,r){N + l,R\)i47h) 

where t^ R are the hopping parameters coupling the left (L) 
and right (R) leads to the array of quantum dots. We are ne- 
glecting the possibility of direct coupling between the left and 
right leads, a simplification of no physical consequences, cor- 
responding to the fact that the array of dots has many of these. 
The approach we are formulating with Eqs. ( |46] | and ( l47b dif- 
fers somewhat from that of Refs. [45.47.1 , but in general terms 
the two approaches are perfectly equivalent. 

The important aspect of the tunneling approach proposed 
in Ref. 1,471 is the need to compute the off-diagonal Green's 
function G]^ {z) for accessing the tunneling properties of 
the system, including in the calculation the coupling to the 
semi-infinite leads. This can be done in many different 
ways ll24ll30ll47ll5Tll52ll53ll . leading in the end to results sim- 
ilar to those obtained using non-equilibrium Green's function 
methods.^**'^^ 

The full Hamiltonian of the problem is the sum of Eqs. ( l22b . 
( l46b . and Among the several ways available to compute 
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G]^ {z), one possibility is to use again the Dyson's equation 
approach. This requires that we know the exact Green's func- 
tion of the leads, before the coupling to the system is estab- 
lished, this is we need to compute the Green's function of the 
problem defined by Eq. ( |46] |. The calculation of the Green's 
function of the leads can be done using the same method we 
used in Sec. Ill Al but now with the boundary conditions 



(lIGilm) = (iV|Gfl|m) =0. 



(48) 



In this case, however, it is much easier to obtain the Green's 
function from the usual definition ([T]). The wave function of 
the electrons is given by (let us consider the left lead only) 



Nl 



- J2 sin[(n-l)0]|n), 



(49) 



where Nl is a normalization length that is taken to infinity 
in the end of the calculation. Using the definition ([T]i and 
the wave function ( |49] l the matrix elements of the resolvent 
is given by the integral in the complex plane (where C defines 
a contour over the unit circle w — e^^) 



1 



.\m+n—2\ 



27ri J (J Plw'^ + zw + Pl 



-dw . 



(50) 



The integral in ( DOb can be evaluated using the same method 
it has been used to solve for the Green's function of a chain 
with periodic boundary conditions'*^, leading to 



Gr 



(51) 

with G™,j l(^) standing for the retarded Green's function, 
'W2 — —X + — and x — E/{2(3l), such that < 
1. The result ( BTI ) is a generalization of the particular results 
given in Ref. [56] [the same is true for Eq. (|52]i1. Repeating 
the same arguments for the right lead we obtain 



2l(3i 



.(l_a;2)-i/2 



w. 



m+n-2N\ 



(52) 

and X = E/{2Pfi). Central to our study are the surface 
Green's functions Gq'§l{E) and G'''j^\i]\i j^i r{E) which we 
obtain from Eqs. ( BTT i and ( |52] |. 

In order to determine G\^ [z) we introduce = Hl + 
Ht + Hr and the full Hamiltonian H = + Vl + Vr, 
the free resolvent is G° = {z - i/")"^ As before, G]^ {z) 
is determined solving the Dyson's equation, which due to the 
short range nature of Vl and Vr has an analytical solution, 
reading 



G,V^(z) = {ll\G^\Nr)D- 



(53) 



with 



D = (l-ii(ll|G°|lOGoo,L) 

X {l-tl{Nr\G''\Nr)GN+iN+i.R) 

- ii4(l/|G°|iVr)(7Vr|G"|lOGoo,LGw+ijv+i,fl(54) 



A formally equivalent result to Eq. ( |53] ) was first derived 
by Caroli et al.,^^ in the context of tunneling across a one- 
dimensional wire, which was the first application of the 
Keldysh^^ formalism to tunneling problems. 

Let us now describe briefly the calculation of the linear con- 
ductance, for which Eq. ( |53T l is needed. The central quantity 
in this approach is the T— matrix. Starting from the Dyson's 
equation (|7]i and introducing an iterative solution, we arrive at 
an equivalent form for G, given by 



G = G° + G°fG° , 



where T is the T— matrix given by 

f = V + VG^f = V + VGV , 



(55) 



(56) 



which describes the scattering of an electron from an initial 
state |i) = \0L) in the left lead, to a final state |/) = |iV + l_R) 
in the right lead. Assuming that the chemical potential dif- 
ference between the to leads (electron reservoirs) is fiR ~ 
jiL + eV, with e the modulus of the electron charge and V the 
electromotive potential between the two reservoirs, the trans- 
mission rate is given by 



- = ^iY.f^^'-l^L)[l-f{E^,-HH) 



(57) 



where 



TaL.N+iR = {mT\N+lR) = {QL\VGV\N+1R), (58) 
leading to a linear current 



3 



a, (3 



df{x) 



dx 



\Tol,n+ir\'^S{E^ — eJ^) . 



(59) 



As usual, the conductance is given hy g = j/V. At low tem- 
perature, g reads 



5 (ml) 



2^ 
27re2 



|7oL,W+li?PPL(ML)Pi?(ML) (60) 
tltl\Gl''{tlL)\'pL{pL)pRO^L)- (61) 



Naturally, the density of states pl{pl) and pr{pl) should be 
interpreted as the local density of states at sites and + 1, 
respectively. Equation ( IMT l is formally equal to that derived 
by Caroli et a/., using the Keldysh formalism, and should be 
multiplied by a factor of two due to the spin of the electrons. 

We should stress here that our approach to the tunneling 
problem is similar to that developed by Mujica et fl/.,iilM2 
but not exactly identical. These authors used the Newns 
coupling**^ of the system to the leads, and their approach to 
the solution of the Green's function is based on the properties 
of tridiagonal determinants after using Lowdin's partitioning 
technique, whereas we explicitly solve the Dyson's equa- 
tion. 
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B. Application to a three-sites quantum dot 

As an application of the formalism that can be worked 
out analytically in full detail, we study the transport across 
a three-sites quantum dot, as that shown in Fig. [3] Although 
this is a very simple model, it is relevant enough for the our 
illustrative purposes. 




Figure 3: A three-sites quantum dot characterized by t/s. This is 
coupled to two metal leads by the couplings Ir and tL- 



particular study of the effect of an adatom on the conductance 
of a quantum wire was done by Kwapiriski."^ The presence of 
the adatom adds an extra term to the Hamiltonian (l22T i of the 
form 

Had = Had + ^ad , (65) 

with 

H% = eQ\AD){AD\, (66a) 
Vad = tQ{\AD)(xl\ + \xl){AD\) (66b) 

where \AD) is the electronic state in the adatom, is the 
electronic hopping between the adatom and the {x,l) site of 
the array, and is the local electronic energy in the adatom. 

It is again straightforward to apply the Dyson equation for- 
malism to compute the exact Green's function in the presence 
of the impurity, leading to 



The Green's functions of the quantum dot depicted in Fig. 
|3]has been computed in Eq. ( |39] |. for the case where the cou- 
pling to the leads was neglected. Next, we need to compute 
Eq. i53[ for this problem. To this end we need the surface 
Green's functions Gqo.l and G22.R, which lead to 

Gl>tL{E)^G^2tRiE) = f^-j:VT^, (62) 

with X = E/{2f3i), with i = L,R. Interesting enough, the lo- 
cal density of states computed from (|62| |. does not diverge at 
the band edge, as it happens with the density of states of an in- 
finite one-dimensional tight-binding model. The calculation 
of the matrix element of the T— matrix leads to 

|G;V'^(ml)P = 16ii/3«(tA +E- e,fV-^ , (63) 
with T) given by 

V = 4(4/32 _ ^2) [(^2^ ^ tl)p\t\ -{E- e^f) 
+t\tlE{E - eo)] ' + [2E[3\tl + t],){-t\ + [E ^ e^f) 
-2tltliE' - 2p'){E ~ eo) 

+ApHtA+E-e„)\2t^-E + e„)]\ (64) 

The algebraic form of V suggests that for an array of quantum 
dots, giving a simple analytical form as Eq. i53[ may not be 
possible in general, and the last steps of a given particular 
calculation may have to be done numerically. 

C. Defects 

In this section we describe the effect of defects on the elec- 
tronic spectrum of the array of dots as well as on its trans- 
port properties. Let us again consider the generic situation 
described in Fig. [T] and consider as a simple and specific ex- 
ample that at the site (x, I) (1 < a; < N) there is an adatom. 
We want to study what is the effect of this adatom on the spec- 
trum of the system and latter on its transport properties. A 



{na\G\mf3) = (na|G°|m/3) + {na\G'^\xl)T{xl\G°\ml3) , 

(67) 

with the T matrix given by 

T ^tl[z~eQ-tl{xl\G°\xl)]-\ (68) 

and the matrix elements {na\G'^\nif3) are computed from the 
resolvent (5° = (z - HTy \ with Ht defined by Eq. (|22|. 
From Equation ( |67] ) we can compute Eq. ( |53] ) and from this 
the conductance given by Eq. dHTT l. 

In addition, we can compute the Green function of the im- 
purity, determining how the energy level is modified due to the 
coupling to the bath of electrons propagating along the array. 
This is given by {Gad = {AD\G\AD)) 

GAD = [z-eQ-tl{xl\G''\xl)rK (69) 

From Equation ( |69l ), the local density of states at the impurity 
can be computed as pad = ^l/i'SG^^|,. The accepted 6 
values are now the solution of 1 /Gad = 0. 

D. Defects: an application 

Again we make a simple application of the formalism of 
Sec. Hire] considering the system depicted in Fig. |4] 

The Green's functions entering in Eq. ( 1531 ) are given in 
this example by 

G^i = {X-l/2){l + 2X){2tAX + Ae)V-\ (70a) 

G°, = [(A- l/2)(l + 2A)(2tAA + Ae) 

- Xtl/tA]V-\ (70b) 

Gl = -{2t/:,X + Ae)V-\ (70c) 

with T) given by 

V = (A-l/2)[2iA(-l + A+2A2)(2tAA+Ae)-<|(l + 2A)] , 

(71) 

with Ae — eo — e^At is interesting to note that the eigenvalue 
A = 1/2 of the non-perturbed triangular dot is not modified 
by the adatom. The calculation of the conductance is now just 
a matter of using Eq. ( iTOl ) in Eq. ( l53T l. 
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Figure 4: A triangular dot with the atoms coupled together by t^, 
which in turn is coupled to two metal leads by the coupling constants 
tR and tL- An add atom (represented by the square) is connected to 
one of the atoms of the dot via the hopping tQ . 



IV. DISCUSSION 

We have presented in full detail a method to compute the 
lattice Green's functions of an array of quantum dots for the 
cases when the array is isolated as well as when it is coupled 
to two metallic leads. The effect of the leads is to produce 
a self-energy which has both a real (SftS) and an imaginary 
(31]) parts. In the case of a single quantum dot, 3fiS will 
renormalize the energy levels in the dot whereas makes 
the energy levels non-stationary. Both terms contribute to the 
charge transport through the dot. 

The formalism is general and flexible enough to allow for 
the study of how localized defects affect both the energy spec- 
trum and the transport properties. We can consider both the 
case when the defect acts as local potential, i.e., diagonal dis- 
order, and when the defect changes locally the values of the 
hopping integrals, i.e., off-diagonal disorder Also more than 
one defect can be attached to the quantum dot array, at differ- 
ent positions in the lattice. In the case of a random distribu- 
tion of impurities an approximate treatment such as the CPA^ 
can be used to compute the full lattice Green's function self- 
consistentlyi2^ 

The generalization of the present approach to two dimen- 
sions should present no difficulties, allowing for the possi- 
bility to proceed analytically in the calculation of the en- 
ergy spectrum and transport properties of finite size two- 
dimensional ribbons. This possibility will be explored in a 
forthcoming publication. The more restrictive aspect of the 
method could be related to the calculation of the allowed val- 
ues of 9, since these are to be computed at the same time the 
values of the energy eigenvalues are determined. For an array 
of N dots, each dot having Ng sites, the determination of the 
spectrum following a brute force approach would require the 
diagonalization of a matrix of dimension {Ng x iV)^. In our 
approach this is reduced to the determination of the zeros of a 
polynomial of degree Ng. If we consider the case of periodic 
boundary conditions, the values are given by 6^ = 2tt£/N, 
with £ ^ 0, ±1,±2, . . .iV/2 (assuming N even), and the 
equation giving the energy spectrum is the same which we 
would have obtained if we had done a Fourier transform in 
the initial Hamiltonian. 

It important to stress that our approach can easily include 
the case where the quantum dot is represented by a continuous 



model. In this case the Green's function of the dot has the 
form G(r, r', E), where r and r' are two-dimensional vectors 
characterizing the position in the dot. In order to apply the 
developed formalism we only need to choose the values of r 
to which the dots connect among themselves. 

Finally we note that our description of the transport is eas- 
ily generalized to include finite values of the potential bias be- 
tween the leads. In this case, an appropriate treatment of the 
problem requires to solve for the Green's function together 
with an iterative solution of the Poisson's equation. For quasi- 
one-dimensional systems this does not require powerful com- 
putational facilities. This will be addressed in a forthcoming 
publication. 
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Appendix A: AN ALTERNATIVE SOLUTION OF THE 
FINITE CHAIN PROBLEM I 

We develop here an approach to the solution of the finite 
chain problem that builds the symmetries between the two co- 
ordinates of the Green function form the start. This method 
can also be used to tackle the more general problem of the 
quantum dot array. In terms of the resolvents and of V, 
Dyson's equation can take two alternative forms, 

G = G" + G°VG, (Ala) 
G = G" + GVG°. (Alb) 

Forming matrix elements with the basis state vectors and us- 
ing i = J2n=i obtains 

{n\G\m) - G"S,,m 

- tG"{{n - l\G\m) + {n + l\G\m)) ,(A2a) 

{n\G\m) = G"Snrn 

- tG°{{n\G\m - 1) + {n\G\m + 1)) j(A2b) 

with G" — (z— eo)^^. By taking the sum and the difference of 
these equations, one derives equivalent conditions which are 
more symmetrical in the two coordinates of the Green func- 
tion, 

■^{n\G\m) - dnrn-l[{n^l\G\m) + {n+l\G\m) 
+ {n\G\m - 1) + {n\G\m + 1)^ , (A3a) 

= {n-l\G\m) + (n+ l|G|m) 

- {n\G\m + 1) - (n\G\m + 1) . (A3b) 
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Apart from the Kronecker delta term 5nm, Eq- ( I A3 al l defines 
the wavefunction of two particles with a tight-binding Hamil- 
tonian and an eigenvalue l/G". It is obvious that the differ- 
ence of two solutions of this set of equations will be a solution 
of the corresponding homogeneous system (without the Snm 
term). Our strategy for finding {n\G\m) = Gnm is the follow- 
ing: (i) we construct the general solution of the homogeneous 
system assuming plane waves in the two "particle" coordi- 
nates; (ii) we find one solution of the full non-homogeneous 
system, by allowing different plane wave solution for n < m 
and n > ni, very much in the spirit of the Bethe solution for 
two particles with local interactions, moving in one dimen- 
sion; (iii) the general solution is just the particular solution 
of (ii) added to a general linear combinations of the solutions 
of (i). To determine the latter we require additional boundary 
conditions, which, for a finite chain are: 



Gom 


= 0, 


(A4a) 




= 0, 


(A4b) 


GnO 


= 0, 


(A4c) 


GnN+1 


= 0. 


(A4d) 



We begin by writing a general solution of the homogeneous 
system in the form 



Ae 



i0in-\-i02m 



(A5) 



Inserting this trial solution in Eq. ( |A3b| ) we get the condition 
cos 6*1 — cos 02 = 0; we must have 9i — ±62 = 0- With this 
condition, t/jnm is a solution of the homogeneous version of 
Eq. d A3 al l provided 



ij--2tcos0. 



(A6) 



So the solution of the homogeneous equations is a linear su- 
perposition of waves 



(A7) 



where 6 solves Eq. iA6\ . 

We now address the determination of one solution of the 
full non-homogeneous Eqs. (|A3l l, which we write in the form 
4>nm = ipnm' for n < TO, and = V'mn, for 71 > m, where 
il^nm ^rid ipnm two different solutions of the homogeneous 
system. There are only two conditions that mix t/j^ and 
namely 



These conditions cannot be fulfilled by solutions which are 
function of n + to, so we must have: 



Inserting these trial functions in Eqs. ( IA8I ) and solving the 
corresponding linear equations for the constants, one gets. 



A-" 



-A> = B' 



1 



4ii sin 61' 



leading to a solution 



1 



■ sm [a \n — to 



2t sm 

exactly as found in section UlAl eq.(fT4]i We have now carried 
out points (i) and (ii) outlined above, and obtained the general 
solution of Eqs. ( IA3I ). as 



Gn 



sin [6 \n — to|) 



2tsin6' 

where is superposition of waves of the form ( IA7b with 
6 satisfying Eq. iA6i . To enforce the boundary conditions 
it proves more convenient to write the solution in sines and 
cosines as 

Gnm — A COS [9 {n — to)] + B sin [6 [n — to)] 
+ C cos [0 {n + to)] + D sin [9 (n + to)] 
sin {6 \n — to]) 
2t sin 6* ' 

To derive the values of these constants, we insert this solution 
in the Eqs. (IA4I) . use the linear independence of the sine and 
cosine functions and arrive at the final result: 



Gnm (^) 



1 COS [0 {N + 1)] 



2t sin [61 (A^ + l)]sin6' 
X {cos [9 {n — to)] — cos [9 {n + to)]} 
1 ( sin [9 {n + m)] sm[9\n — m\ 
2t 



A9) 



which is the same solution as Eq. ( [TtI i. 



Appendix B: AN ALTERNATIVE SOLUTION OF THE 
FINITE CHAIN PROBLEM II 



Q^i^nn = 1-2 (^n+ln + "^n-ln + V4„+l + V',m-l) > 
^^n+lTi = — 2 i'^-'n+2n + V'nn + ^n+ln+1 + ''Pn+ln-l) ■ 

Because tp'^ and tp^ are solutions of the homogeneous sys- 
tem, we can easily transform these conditions into 

1 = ^ {'>Pn+ln - ^n+ln + 4'nn~-l ^ "^rm-l) , (A8a) 
= V-™ - i'nn + V-^+ln+l ^ V^n+ln+l • (A8b) 



In the previous appendix, the Dyson equation was written 
in two alternative forms, see Eqs. (lAll) and (IA2I) . 

In Eq. ( IA2ab . the ket is unchanged and it can be thought 
of as a tight-binding equation for the bra of {n\G\m) with an 
inhomogeneity at site to. Since we are dealing with a real 
Hamiltonian, we can make the following ansatz for Gn m — 
(nlGlTO): 

Q _ / Gnm — '^1 ("^) COS 9n + 02 (to) sin 9n , n < m 

(Bl) 
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where ai{m) and bi{m) (i = 1,2) are arbitrary functions of where Ci{n) and di{n) (i ^ 1,2) are arbitrary functions of n. 
m. 

In Eq. ( IA2bl i. the bra is unchanged and it can be thought of 
as a tight-binding equation for the ket of {n\G\m) with an in- 
homogeneity at site n. We can thus make the following ansatz: 



Gnm — ci [n) COS 9m + C2 {n) sin ( 



di (n) cos 9m + d2 (n) sin 9m 



m > n 
m < n 
(B2) 



Combining Eqs. dBlb and ( IB2I ). we arrive to the following 
ansatz for the Green's function: 



J 



Gri,7n, — 



fli cos 9n cos 9m + a2 cos 9n sin 9m + aj, sin 9n cos 9m + 04 sin ^^r^ sin 9m , n < m 
= bi cos 9n cos (?m + &2 cos 0n sin (?m + 63 sin 0ri cos 9m + 64 sin 0n sin 9m , n > m 



nm 



r 



(B3) 



Notice that now the coefficients are site-independent. 

From the boundary conditions G^„j = G^q = 0, we obtain 
fli = 02 = 61 = 63 = 0. From the boundary conditions 
^N+irn = = 0' ^6 have = -04 tan 6I(7V + 1) and 

&2 = — ^4 tan 9{N +1). The matching condition atn — m — 
1 yields = (continuity of the Green's function) 

and thus 04 = 64. The matching condition ain = m yields 
^m-im ~ ^m-im — (discontinuity of the derivative of 
the Green's function) and thus 04 — {sm9 tan9{N + l)i)~^. 
The final result can therefore be written as 



with the functions Li, with z = 0, . . . , 5 given by 



G 



< 

nm 



1 — tan 0(iV + 1) sin 9n cos 9m + sin 9n sin 6*™ 
7 



tan6'(iV + l)sin( 
1 — sin 9n 



1 sm9{N 
1 sin0(7V + l)sin6' 



G 



> 

nm 



G 



< 

mn 



(B4) 
(B5) 



Again 9 is determined by the dispersion relation Eq. ( IA6I 1. 
This yields an alternative (but equivalent) representation of the 
Green's functions of a tight-binding chain with open bound- 
aries. 

The extension to the more general case is analogous, but 
one has to take special care by defining the matching condi- 
tions because the unperturbed Green's function is now a ma- 
trix. It then follows that the Green's function for the non- 
diagonal matrix elements which are not constrained by the 
boundary conditions will be discontinuous for energies which 
are not eigenenergies of the unperturbed system. 



Appendix C: MATRIX M OF EQ. OH) AND MATRIX V OF 
EQ.<^ 



Lo 


= cos{m9) 


- cot(A^6l + 9) sin{m9) , 


(C2a) 


Li 


= cos{m9 - 


- 9) - cot(iV6' + 9) sm{m9 H 


- 0) (C2b) 


L2 


= [cot {m9) 


- cot(7V6' + 9)] sm{m9 - 9 


), (C2c) 


L3 


— siii{m9 - 


-0), 


(C2d) 


Li 


= cos{m9 - 


- 9) - cot{N9 + 9) sm{m9 H 


^ 0) j(C2e) 


^5 


= sin(m0) . 




(C2f) 



The matrix V of Eq. ( [33] l is given by 



V 



cos{m9) sm{m9) 
cos(m6'-6l) sin(TO6l - 6*) 
cos(m6i + 9) a{9) b{9) 



(C3) 



with the functions a{9) and b{9) given by 

a{9) = cos{m9 + 9)tan{m9) , (C4a) 
b{9) = sm{m9 + 9)- cos{m9 + 9) tan(m6') . (C4b) 

The vector q-'^ = [qi, 92, 93] entering in Eq. (l33b has its com- 
ponents given by 

51 = {tPL)-^[tGliPL~t^GliGl^Glsin{m9 -9) 
X sin{N9 + 9 -m9) 

+ t^G°„Glsm{N9~m9)P.,r,-i)], (C5a) 

52 = {tPLy^[-~PL + t^Gl^Gls\xv{N9 + 9-m9) 

X s\TL{m9-9), (C5b) 

93 = {tPLy^[-PL + t^Gl^Gls\xv{N9-m9) 

X sin(TO6'). (C5c) 



The matrix M of Eq. ( [30l ) is given by 



M = 



Lo + tG^^Li tG^L2 

tG^j.L4 L^ -f- (,i_r[j,iJ3 

tG'^^.Li tG'i^{L2 — L4) —tGfiL^ 





tG^^L^ 



(CI) 



Appendix D: FULL ANALYTICAL EXPRESSIONS FOR THE 
GREEN'S FUNCTIONS 



After using the boundary conditions and three of the four 
time reversal conditions, the Ansatz for the Green's functions 
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for n < m is 



Gnm 
U 



Gnm 
r' 



— [cos{n9) — cot{m9) sm{nd)] 
+ A> [cot(m6i) - cot(7V6' + 9)] sm{n9) , (Dla) 
= A>[cot{m9) - cot{N9 + 9)] sm{n9) , (Dlb) 
B<sm{n9), (Die) 
(Did) 



G»™ = ^<cos(n6 



sin(n6 



and for n > mis 



Gnm 

II — 

Gnm 

rl — 

Gnm 

Ir — 

r^nm 

rr 



A>[cos(n9) - cot{N9 + 9) sin(?i6i)] , 



A>i cos{n9) 



B> sm{n9) , 



A^^[cos{n9) - cot{L9 + 9) sm{n6 
A>^[cos{n9) - cot{ni9) sin(n6')] 
B^j. sin{n9) . 



(D2a) 
(D2b) 

(D2c) 



Following the method described in the bulk of the paper, the 
full analytical expressions for the Green's functions are given 
by 



and 



sin{i9) sin[(iV + 1 - j) 



i < j 



(D5) 
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s\-a{j9) sin[(7V + 1 - i) 



i> 3 



(D6) 



Note that Gj^ = G^; and that the diagonal Green's functions 
obey 
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(D8) 



Qij 



^11 



-^sin{i9)PL_j, i<j, 

G 

-^sm{j9)PL-t, i>j, 



-^sin[(iV + l-j)^]P,_i, z<j 
^Gj^sin[{N+l-i)9]P^^^, i>j 



(D3) 



(D4) 



which is similar to Eqs. (|13b| i and (|32a| i. which is the general- 
ization to the lattice of the discontinuity of the first derivative 
of a Green's function. 

Similar results to those given in this Appendix have been 
also obtained in Ref. |58] in the context of organic molecular 
systems, but no hints about the method used to derive them 
was given. 



* Kenichi Nishi, Device Applications of Quantum Dots, chap. 12. 
In Y. Masumoto, T. Takagahara, editors. Semiconductor Quan- 
tum Dots: Physics, Spectroscopy and Applications, (New York: 
Springer, 2002) p. 457. 

^ Charles Z. Hotz and Marcel Bruchez, editors. Quantum Dots: Ap- 
plications in Biology , (Totowa: Humana Press, 2007). 

' Michael C. Petty, Molecular Electronics: From Principles to 
Practice, (Ames: WileyBlackwell, 2007). 

" T. Kwapinski, J. Phys.:Condens. Matter 18, 7313 (2006); idem, 
ibidem 19, 176218 (2007). 

^ Wenjie Liang, Marc Bockrath, Dolores Bozovic, Jason H. Hafner, 
M. Tinkham, Hongkun Park, Nature 411, 665 (2001); Jing Kong, 
Erhan Yenilmez, Thomas W. Tombler, Woong Kim, Hongjie Dai, 
Robert B. Laughlin, Lei Liu, C. S. Jayanthi, and S. Y. Wu Phys. 
Rev. Lett. 87, 106801 (2001). 

* Alexander Onipko, Yuriy Klymenko, and Lyuba Malysheva, J. 
Chem. Phys. 107, 5032 (1997). 

^ Alexander Onipko, Yuriy Klymenko, and Lyuba Malysheva, Phys. 

Rev. B 62, 10480 (2000). 
^ Mario I. Molina, Phys. Rev. B 67, 054202 (2003); idem, ibidem 

73, 014204 (2006). 
' idem, ibidem 74, 045412 (2006). 

D.A. Ryndyk, R. Gutierrez, B. Song, and G. Cuniberti, 

arXiv:0805.0628vl. 
" Zhao Yang Zeng, Yi-You Nie, F. Claro, and W. Yan, Phys. Lett. A 

331, 84 (2004). 

Hua Li, Tianquan Lu, and Punan Sun, Phys. Stat. Sol. (b) 242, 
1679 (2005). 



" Liliana Arrachea, Phys. Rev. B 72, 125349 (2005). 
Akira Oguri, Phys. Rev. B 63, 115305 (2001). 
C. Karrasch, T. Enss, and V. Meden, Phys. Rev. B 73, 235337 
(2006). 

Hiroyuki Ishii and Takashi Nakayama, Phys. Rev. B 73, 2353 1 1 
(2006). 

Frank Reininghaus, Thomas Korb, and Herbert Schoeller, Phys. 
Rev. Lett. 97, 026803 (2006). 

C. Caroli, R. Combescot, P. Nozieres, and D. Saint-James, J. Phys. 

C:Solid St. Phys., 4, 916 (1971); idem, ibidem, 4, 2598 (1971); 

idem, ibidem, 5, 21 (1972). 
" R. Combescot, J. Phys. C:Sohd St. Phys., 4, 2611 (1971). 

R. Combescot and G. Schreder, J. Phys. C:Solid St. Phys., 6, 1363 

(1973); idem, ibidem, 7, 1318 (1974). 

F Guinea and J. A. Verges, Phys. Rev. B 35, 979 (1987). 

P Sautet and C. Joachim, Phys. Rev. B 38, 12238 (1988). 
" H. Mizes and E. Conwell, Phys. Rev. B 44, 3963 (1991). 

N. M. R. Peres and F Sols, J. Phys.: Condens. Matter 20, 255207 

(2008). 

A. R. Akhmerov, J. H. Bardarson, A. Rycerz, and C. W. J. 
Beenakker, Phys. Rev. B 77, 205416, 2008. 
N. M. R. Peres, F. Guinea, and A. H. Castro Neto, Phys. Rev. B 
73, 125411 (2006). 
" N. M. R. Peres, A. H. Castro Neto, and F. Guinea, Phys. Rev. B 
73, 195411 (2006). 

A. H. Castro Neto, F. Guinea, N. M. R. Peres, Physics World, 
November, 33 (2006). 

Katsunori Wakabayashi, J. Phys. Soc. Jpn. 71, 2500 (2002). 



13 



B. H. Teng, H. K. Sy, Z. C. Wang, Y. Q. Sun, and H. C. Yang, 
Phys. Rev. B 75, 012105 (2007). 
" David K. Ferry and Stephen M. Goodnick, Transport in Nanos- 
tructures, (Cambridge: Cambridge University Press, 2001). 
Supriyo Datta, Electronic Transport in Mesoscopic Systems, 
(Cambridge: Cambridge University Press, 2005). 
Eleftherios N. Economou, Green's Functions in Quantum 
Physics, (Berlin: Springer- Verlag, 2006). 
G. Barton, Elements of Green's Functions and Propagation: Po- 
tentials, Diffusion and Waves, (Oxford: Oxford University Press, 
1989). 

Shigetoshi Katsura, Tohru Morita, Sakari Inawashiro, Tsuyoshi 
Horiguchi, and Yoshihiko Abe, J. Math. Phys. 12, 892 (1971). 
Katja Lakatos-Lindenberg, Richard P. Hemenger, and Robert M. 
Pearlstein, J. Chem. Phys. 56, 4852 (1972). 
" Peter Karadov and Obis Castano, J. Mol. Struct. 90, 115 (1982); 
idem, J. Chem. Soc. Faraday Trans. 2, 78, 73 (1982). 
Robert Vein and Paul Dale, Determinants and Their Applications 
in Mathematical Physics , (New York: Springer- Verlag, 1998). 
K. S. Dy, Shi-Yu Wu, and T. Spratlin, Phys. Rev. B 20, 4237 
(1979). 

A. A. Bahurmuz and R D. Loly, J. Math. Phys. 22, 564 (1981). 
WiUiam A. Schwalm and Mizuho K. Schwalm, Phys. Rev. B 37, 
9524 (1988). 

'^'^ Milton Abramowitz and I. A. Stegun, Handbook of Mathematical 
Functions, (New York: Dover, 1965). 

'^^ Ronald Bass, J. Math. Phys. 26, 3068 (1985). 

Rolando Perez-Alvarez and Federico Garcia-Moliner, Transfer 
Matrix, Green Function and Related Techniques, (Castello de la 
Plana: Publicaciones de la Universitat Jaume I, 2004). 



D. M. Newns, Phys. Rev. 178, 1123 (1969). 

^ Sydney G. Davison and Kenneth W. Sulston, Green Function The- 
ory of Chemisorption, (Dordrecht: Springer, 2006) chap. 4. 

^'^ V. Mujica, M. Kemp, and M. A. Ratner, J. Chem. Phys. 101, 6849 
(1994); idem, ibidem 101, 6856 (1994). 

'^^ Vladimiro Mujica, Mathieu Kemp, Adrian Roitberg, and Mark 
Ratner, J. Chem. Phys. 104, 7296 (1996). 

Vladimiro Mujica, Mathieu Kemp, Adrian Roitberg, and Mark 
Ratner, A model for Coulomb interaction in electron transport in 
one-dimensional mesoscopic devices. In E. Ludena, P. Vashishta, 
R. Bishop, editors. Condensed Matter Theories: Volume 11, 
(Hauppauge, New York: Nova Sciences Publishers, 1996) pp. 
261-272. 

'° Daniel S. Fisher and Patrick A. Lee, Phys. Rev. B 23, 6851 (1981). 
" T. Ando, Phys. Rev. B 44, 8017 (1991). 

P. A. Khomyakov, G. Brocks, V. Karpan, M. Zwierzycki, and P. J. 
Kelly, Phys. Rev. B 72, 35450 (2005). 



K. Nikolic and A. MacKinnon, Phys. Rev. B 50, 11008 (1994). 



53 

Antti-Pekka Jauho, Ned S. Wingreen, and Yigal Meir, Phys. Rev. 
B 50, 5528 (1994). 
'5 J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 

Z. G. Yu, D. L. Smith, A. Saxena, and A. R. Bishop, Phys. Rev. B 
59, 16001 (1999). 

Per-Olov Lowdin, J. Math. Phys. 3, 969 (1962). 
Alexander Onipko and Lyuba Malysheva, Coherent Electron 
Transport in Molecular Contacts: a Case of Tractable Modeling, 
chap. 23. In Sergey Edward Lyshevski, editor, Nano and Molecu- 
lar Electronics Handbook, (London: CRC, 2007). 



